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Abstract: A strategy for calculating the color structure needed for soft gluon resum- 
mation for processes with any number of colored partons is introduced using a iV c — ► oo 
inspired basis. In this basis a general formalism can be found at the same time as the 
calculations are simplified. 

The advantages are illustrated by recalculating the soft anomalous dimension matrix for 
the processes gg -> gg, qq -> qqg and qq -> ggg. 
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1. Introduction 

The strong force comes with the problem of being precisely strong. Although the coupling 
constant of QCD is small enough for perturbation theory to make sense at all, it is large 
enough to cry out for higher order corrections for many processes, and in some regions of 
phase space, large enough to invalidate a fixed order calculation. 

This is the case in the collinear region, where a large logarithm compensates for the 
moderate smallness of a s , and similarly in the soft region where there is a large effective 
phase space ~ log(hard scale/soft resolution scale) in transverse momentum. In these re- 
gions resummation methods are needed. In the collinear DGLAP region [1—4], where the 
emission can be seen as coming from one parton, the color structure is trivial and Sudakov 
form factors can be used to describe no-emission probabilities. 
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Unfortunately, the strong force is not only strong, it is also complicated, in the sense of 
being non-Abelian. In the soft region, where emissions have contributions from branchings 
of different partons this complicates matters. The real emission coming from the inter- 
ference term of emission off parton i and emission off parton j, is canceled by the virtual 
gluon exchange between parton i and j. (Using Feynman gauge self-energy type diagrams 
can be neglected, i.e. i ^ j.) 

Under the assumption that emissions strongly ordered in transverse momentum domi- 
nate, all leading logarithms in (hard scale/soft scale) from virtual corrections exponentiate 
and can be resummed. However, since these gluon exchanges affect the color structure, the 
exponentiation must be done at the amplitude level. Thus a no-emission amplitude 

I Q \ 

-ll a ° {k ^ )d 4t T \ m ^ (L1) 

\ Qo J 



M = exp 



can be derived. In the above, Mo is the undressed hard scattering amplitude as a vector 
in color space and T is a matrix in color space, describing the effect of exchanging gluons 
between the various partons, 

r = J]%c«. (1.2) 

i<j 

Here C 13 describes the color algebra part and fi,-,- contain the azimuth and rapidity mo- 
mentum integral over the exchanged gluon k', 



n = _l ( _ 1) l f dy'dfi k '±Pi-Pj j (1 
13 2 { ' J n 2tt 2 Pi -k'k'- Pj 2 { 



(1.3) 

■Pj " 

with Sij = — 1 if the partons ij are both incoming or both outgoing, and 1 otherwise, and I 
counts how many of the involved partons which are quarks in the initial state, anti-quarks 
in the final state or gluons, assuming the convention in Eq. ([O]) for the triple gluon vertex. 
In the above equation the «7r-terms, coming from Coulomb gluon exchange, would give rise 
to an unobservable phase in an Abelian theory. For a non-Abelian theory they do, however, 
enter in a physically relevant way. 

In general the color basis used need neither be orthogonal or normalized. In fact, it 
will be seen below that the calculations simplify significantly in a special basis which is 
not. For a non-orthonormal basis, the matrix of scalar products S, calculated by summing 
over quark, anti-quark and gluon indices a,b,c, ... 

s mn =< c m ,c n >= C2 c ..{C n Y abc .., (1-4) 

a,b,c... 



is needed. (Note that C m above is a basis tensor in color space, whereas Q l] in Eq. fll.2p are 
matrices in this basis, describing the effect of gluon exchange between parton i and j.) The 
physical no-emission probability is given by a = JVFSM. As an aside it is pointed out that 
scalar products between tensors corresponding to linear combinations of color structures 
of Feynman diagrams with real coefficients are real. 
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In the simple case of q\q2 — > (73^4 the vector space containing the color structure has 
only two dimensions, and the basis vectors are often taken to be the "i-channel singlet 
octet basis", 



^91929394 ~~ ^9194^9293 
f 2 — t 9 t 9 — - 

91929394 9491 9392 2 



"9193^92 94 ^9194^9293 



J1.5) 



In this case the issue of keeping track of the color structure amounts to a moderate com- 
plication. However, already for gg —> gg & six dimensional vector space is needed (reducing 
to a five dimensional space for N c = 3) and for gg — > <?<7<? there are 22 different color states 
to keep track of (reducing to 16 for N c = 3) [5,6]. In the later case, to keep track of the 
change in color structure as a result of virtual gluon exchange between a pair of partons, 
one naively - without using further symmetries, thus needs to calculate the effect of gluon 
exchange on 22 different color states, and then decompose the result into the 22 different 
color tensors by taking scalar products, implying in total 22 2 scalar product. (This number 
may be reduced, for example by using the fact that the soft anomalous dimension matrices 
are symmetric if stated in orthonormal bases [7,8].) The color structure thus gives rise to a 
major computational complication, and so far the soft anomalous dimension matrices have 
only been calculated for the 2 — > 2 processes [5,9-14] and the 2 — ► 3 processes [6, 15]. (For 
observable related and experimental work, see for example [16-24].) 

If one is only interested in a fixed order expansion, for example as in [25], there is 
no need to choose an explicit basis. Indeed the soft anomalous dimension matrix can be 
written down in a compact basis-independent way for any number of partons, both at 
one-lopp and two-loop order [42]. Similarly, for the purpose of deriving general theoretical 
properties it is often wiser to stay basis independent, and several interesting results have 
recently been derived without explicit basis choices [8,26-30]. However, to actually perform 
the numerical exponentiation of Eq. (|1.1|), to obtain all-order results, an explicit basis is 
needed. 

It is thus clearly desirable to find a simplifying general strategy. Especially, a unified 
formalism is needed for the long term goal to incorporate non-leading color effects in event 
generators. The major current event generators all work in the leading iV c limit [31-35]. 
This means that the color structure is decomposed into leading N c contributions, using 
Eq. (|2.2| ) and Eq. (2.1) below. Color suppressed interference terms between different color 



structures are neglected. It was argued a long time ago that for gluon amplitudes with 
fixed power of a s these terms are suppressed by 1/-/V 2 [36]. However, there may in general 
be many suppressed terms. As an example consider N g — 2 gluons attached in a row to 
one gluon line, giving in total N g gluons. (For N g up to five, all tree level graphs have this 
topology.) The squared amplitude is given by 

Nc 9 ~ 2 {Nl - 1). (1.6) 

If the diagram is decomposed into different color topologies (which are orthogonal in the 
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N c — > oo limit) the sum of the parts squared separately is 



[(^-i^ + ^.^-l)]. (1.7) 



When N c — * oo both expressions grows as A^ 9 and their ratio approaches one. However 



for finite N c the difference grows with N g and already for N g = 4, if N c = 3, Eq. Ql,7| ) is 
only 19/27 of Eq. (O), [37]. For 7 gluons Eq. (O) is less than 50% of Eq. (O). 



The method suggested in this paper for dealing with the color structure of multi- 
parton processes is developed with the resummation of soft gluons in mind, but clearly, as 
it describes the effect of gluon exchange on any colored amplitude, it may also prove useful 
for NLO (and higher order) corrections to amplitudes with (many) colored partons. 

The results may also be used to calculate effects stemming from the non-global nature 
of most observables, the ordinary "non-global logs" [38,39], as well as the color suppressed 
"super leading logarithms" carrying extra powers of log(hard scale/soft scale), suggested to 
enter at order af in perturbation theory [25,40,41]. Indeed, as the non-global logarithms 
originate from real radiation outside an experimental exclusion region, to calculate the 
contribution from n emissions outside the exclusion region requires the soft anomalous 
dimension matrices for processes containing n additional partons. 

As the two-loop soft anomalous dimension matrices have been proven to be propor- 
tional to the one-loop results (for processes with any number of colored and uncolored 
massless external legs), the present method can trivially be used also for two- loop anoma- 
lous dimension matrices [42]. Recently it has been suggested that similar results also hold 
for the three- loop anomalous dimension matrices and that they may hold to any order, as 
long as the partons remain massless [26-28]. For massive external legs this simple relation 
breaks down [29,30]. 

The layout of this paper is as follows: First the formalism for constructing a basis is 
described in section |2| and computational rules for gluon exchange in this basis are derived 
in section |3[ To illustrate the advantages with the constructed bases, the soft anomalous 
dimension matrices for gg — > gg, qq — > qqg and qq — > ggg are recalculated in section |||. 
Finally some concluding remarks are made in section |5|. 



2. General basis formalism 

2.1 Construction of a general basis 

Previous strategies for dealing with the color structure needed for resummation of soft 
gluons have lately been based on multiplet decomposition for finding a basis [6, 14]. In 
this way symmetry properties are exploited to construct complete orthogonal bases (which 
easily can be normalized). Clearly, using an orthogonal basis has advantages. The result 
is easy to interpret and the matrix of scalar products between basis vectors is diagonal. 

A disadvantage of the multiplet strategy is, however, that increasingly complicated 
projection operators need to be used, and no closed form for deriving these projection 
operators exists (to the knowledge of the author). 
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Another complication is that the projection operators, which tend to be expressed in 
terms of the symmetric and anti-symmetric structure constants f abc and d a bc, need increas- 
ingly complicated computational rules for contraction of indices, that is, computational 
rules involving more and more /'s and cf s. Alternatively, the structure constants can be 
reexpressed in terms of the generators of the fundamental representation tq iq2 , 




= 2(Tr[rt c ] T Tr[t b tH c }) = 2{t a qiq / q2q ^ qi T CWU). ^ 



In this case any scalar product, of arbitrarely complicated color tensors, can be calculated 
using the gluon index contraction relation 

t 9 cat 9 db = \{5 ad hc-^ ac 5 bd ). (2.2) 

However, the expression for the color structure tensor will contain 2^ °f / s an< ^ ^ s ) terms, 
and the scalar product of the tensor with itself thus 2 2 ^ f s anc ^ ^ s ) terms which each 
has to be contracted separately. 

An alternative strategy would be to construct a basis by starting from a sufficient 
number of arbitrarily chosen color tensors, or by exploiting possible symmetries. This will 
work well for a small vector space, cf. [15], but will tend to give very lengthy expressions 
for the basis vectors if Gram-Schmidt orthogonalization is used for a large vector space. 
On the other hand, if the basis vectors are not made orthogonal the decomposition of color 
structures resulting after gluon exchange will in general be cumbersome. (This complication 
is circumvented in the special non-orthogonal basis suggested below). In addition it has to 
be proved that the basis actually span the relevant space. 

These issues make it worth exploring other strategies for constructing the basis in the 
general case of any number of colored and uncolored partons. The basis clearly has to span 
the relevant space. It may seem desirable to find an orthogonal (normalized) basis, but it 
will be seen below that using a special non-orthogonal, non-normalized basis significantly 
diminishes the computational effort, mainly since the state obtained after gluon exchange 
is immediately, i.e. without taking scalar products, a linear combination of basis states. 
There is thus no need for calculating A^ asig 2 scalar products for every possible gluon 
exchange. 

The solution is to use a basis inspired by the N c — > oo limit. In the case of infinitely 
many colors, two color lines in a Feynman diagram are never the same, and gluons may be 
represented by two color lines going in opposite directions. In this case, all possible color 
structures can be represented by all ways of connecting incoming and outgoing color lines. 
The strategy suggested here is thus similar to methods used in [43-46]. Especially it is noted 
that the bases suggested here for resummation are similar to the color structure treatment 
suggested in [46] to deal with real parton emission in event generators. For N c = oo 
the scalar product between different color topologies, divided by the scalar product of a 
topology with itself, equals zero. However, for finite N c there are scalar product terms 
which are suppressed only by 1/N C . 
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Another important property of the bases constructed in the aforementioned way is that 
they are completely democratic w.r.t. different quarks, different anti-quarks and different 
gluons. This implies, for example, that once the effect of gluon exchange between the 
gluons gi and 32 has been calculated, the effect of gluon exchange between any other 
gluons can be obtained by relabeling of indices i.e. renumbering of basis tensors. One 
therefore never needs to calculate more than six different exchanges gg, qq, qq, gq, gq and 
qq. In addition, it will be seen below that the color structure after a gluon exchange on 
a given color topology is a linear combination of at most four different basis tensors. The 
soft anomalous dimension matrices will thus be relatively sparse in the suggested bases, 
which should simplify numerical exponentiation. 

It is also worth stressing that the suggested bases are well suited for comparison to 
the N c — > 00 limit, as the bases are easy to interpret and the soft anomalous dimension 
matrices will turn out to be diagonal in this limit. This implies that they are ideal for 
comparison to the radiation pattern obtained from event generators tending to work in the 
N c -» 00 limit [31-35]. 

The reduction in calculational effort for the soft anomalous dimension matrix with the 
suggested basis is thus threefold. There is no need to calculate scalar products, reducing 
the computational effort with a factor ~ ^ from the number of scalar products and a 
factor 2^ °f / s an< ^ s ) from the number of terms in each of the scalar products, assuming 
Eq. ( |2.1| ) is used. Furthermore there are at most six, as compared to Np(Np — l)/2 for 
Nr> external particles, different gluon exchanges to keep track of, the others are related by 
relabeling of indices. 

Unfortunately this does not quite remove the bad scaling of the problem with the 
number of partons, as instead of having to calculate ~ -^basis sca ^ ar products for each 
contribution to the soft anomalous dimension matrix, one has to calculate ~ -^basis sca ^ ar 
products between the basis vectors, as they are only orthogonal in the N c — > infinity limit. 
However, this only has to be done once. In addition calculating scalar products using 
Eq. pj| ) and Eq. ( |2~^ ) gives just one, as opposed to 2^ °f / s an d d s ) ) different terms. 

What remains is thus a scaling of type basis' Very roughly speaking iV^^jg tends to 
grow as -/Vp!, cf. section 2.2-2.4. But, bearing in mind that only the topology of the color 
contraction, and not the labeling of indices is important for the scalar product, should 
naively reduce the (iVp!) 2 scaling by a factor ~ N g \N q \Nq\ from the number of ways of 
labeling the indices. What remains is then a factorial growth for processes with only 
gluons. 

Note however, that for processes with many enough external partons, the major com- 
putational effort will not lie in finding an expression for the soft anomalous dimension 
matrix, but in numerical exponentiation of the obtained result. As numerical matrix ex- 
ponentiation scales with the cube of the matrix size, and the number of basis vectors tends 
to grow factorially with the number of partons, calculations with more than ten particles 
seem unlikely. For practical implementations, it is also worth pointing out that the num- 
ber of basis vectors highly depend on the kinds of partons involved. For processes with 
no external gluons and N q = Nq = Np/2 partons, the number of basis vectors is (iVp/2)! 



-6- 



whereas for processes with only external gluons the size of the basis tens to grow rather as 
Npl/e, cf. section |J and [D|. 

That a basis constructed in the above described way is complete for N c = oo is clear 
from the fact that it represents all possible color topologies. For finite N c , some of the 
color tensors may be linearly dependent, and the basis over-complete, but it will still span 
the space. One way of thinking of the reduction in dimension of the color space is to note 
that tensors corresponding to multiplets which are anti-symmetric in more than N c quark 
indices are not possible. Requiring that a color decomposition should be valid for all N c 
defines a unique decomposition of a N c = 3 tensor. 

Another way of convincing oneself that the above bases are complete, is to note that 
every internal gluon line in any Feynman diagram can be removed by first using Eq. ( |2.1| ) 



to remove the triple gluon vertices and then Eq. (2.2) to remove gluon propagators. In 
this way any Feynman diagram, tree level or not, can be decomposed into color structures 
containing no gluon propagators. What remains is a linear combination of color structures 
containing internal quark lines, external quarks, external anti-quarks and external gluons. 
That is, a linear combination of terms of precisely the form obtained by first splitting all 
gluons to qq-paris, and then connecting quark and anti-quark lines in all possible ways. 

Below, the construction of basis tensors will be investigated in more detail, first in the 
special case of external quarks only, then for external gluons only, and finally in the general 
case of both. 

Before moving on we note that from the color algebra point of view there is no difference 
between an outgoing quark and an incoming anti-quark, from here on simply collectively 
referred to as quark, or an incoming quark and outgoing anti-quark, from now on referred 
to as anti-quark. Opposite conventions may be used elsewhere. In addition the placing of 
quark and anti-quark indices on the fundamental generators may be varied. 

2.2 The quarks only case 

Finding a basis in the case of only external quarks is trivial. The basis just consists of all 
possible ways of connecting quarks and anti-quarks. For N q = Ng quarks (clearly, for each 
incoming quark line there is also an outgoing) this can be done in 

%asis = N q \ (2.3) 

ways. The squared norm of these basis vectors, calculated using Eq. ([TJ), is equal to iV c 9 . 
To denote the tensors the notation 



(<?1<?3) (<? 2 <?4) = <W<W (2-4) 

is used. A complete basis for q±q2 — » 9354 is thus the tensors («7i<3 , 3)(<72 < 74) an d ( Q1Q4) (Q2Q3) ■ 
In fact this is the basis used in [10]. 

2.3 The gluons only case 

To construct the basis in the case of gluons only, closed quark loops with external gluons 
attached are used. For example, for four gluons, all gluons may be connected to the same 
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quark line giving (4 — 1)! = 6 topologically different diagrams. Alternatively the gluons 
may be connected two and two in three different ways. Indeed the color space also has 
nine dimensions, however, only half of the linear combinations of the six fully connected 
topologies are physical, due to the fact that quarks and anti-quarks enter QCD on equal 
footing. Therefore, if, in a quark loop, a quark is going around in one direction, the 
topology with the quark going around in the opposite direction (i.e. the gluon index order 
is reversed) must also contribute. 

More explicitly, introducing the notation 

(9l92...9N g ) = Tr[t^...t^] = « 2 V"C g % 1; (2-5) 

to denote N g gluons attached clockwise in the order g\...gN g on a quark line, we note that 
the physical linear combinations must be 

(gm-9N g ) + (-l) Ns (gN g -929i). (2.6) 
To understand the sign, decompose any tree level Feynman diagram with only gluons using 



Eq. (2.1) and Eq. ( |2.2| ). The result is a sum of color structures where the N g gluons are 
attached in different orders to the quark-line. For a specific order, the anti-cyclic order is 
obtained by reversing the direction of the quark-line in every vertex, i.e. taking the other 
term in Eq. ( |2.1| ) everywhere. This gives a factor (— l)^ -2 as there are N g — 2 vertices, 



explaining the sign in Eq. (2.6). 

Thus, in the case of gg — > gg, only six color tensors are needed (for general N c ). This 
explains the observation that some tensors decouple for gg — > gg and gg — > ggg [5,6, 12]. 

The problem of constructing the -/V ff -gluon basis in the general case thus boils down 

to: 

(1) Find all the ways of grouping the N g gluons such that each group contains at least 
two gluons. (Groups with only one gluon would correspond to the color structure 
t 9 qq = 0.) For four external gluons the possible groupings are thus {4} and {2,2}. 

(2) For each fully connected grouping, such as {4}, find all physical different ways of 
arranging the gluons. For N g gluons this gives (N g — l)!/2 different color tensors 
where the factor 1/2 is present since only one combination of the cyclic and anti- 
cyclic ring is physical. 

(3) For disconnected groupings, such as {2,2}, 

(3a) Find separately, for each subgroup, all physically different ways of arranging 
the gluons. 

(3b) Distribute the gluon indices {gi----gN g } in all possible ways among the different 
subgroups. 

(3c) Combine the different sub-groupings in all possible ways, taking into account 
that, if all gluon indices are equal, two groupings do actually correspond to 
the same physical state. For example the subgrouping {{gi, 52}, {93, 9i}} arid 
{{53,54}, {91,92}} are equal. 
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Following this recipe a complete basis describing the color structure for any number of 
external gluons can be constructed. 

Neglecting the issue of physical linear combinations, the possible color tensors coincide 
with the color tensors obtained by replacing each gluon with one quark and one anti-quark 
line, with the important exception that contractions between a qq pair corresponding to 
the same gluon are disallowed. The problem of finding all such topologies is equivalent to 
the number of ways of mapping N elements to each other without mapping a single one to 
itself, which has a known solution 

wffcl)!^. (2 . 7) 

The convergence to N\/e is very quick, rounding off to the closest integer works already 
for N = 1. 

Note however, that this just gives the total number of linearly independent color ten- 
sors (for N c = oo). As mentioned above, only tensor combinations where quarks and 
anti-quarks enter on equal footing are physical. For every quark ring participating in 
building up a color tensor, the corresponding anti-quark ring has to be added. This re- 
duces the number of physical tensors of a certain topology, such as {3, 2} with a factor 
(l/2)# rin § s buildin S U P the tenSOT , that is (1/2) 2 for {3,2}. 

As the number of fully connected color topologies, where all N g gluons are attached to 
the same quark- line equals (N g —1)\, the fraction of color tensors corresponding to fully con- 
nected diagrams is roughly e/N g , again ignoring the issue of physical tensor combinations. 
Tree level QCD Feynman diagrams with only external gluons (more generally, no gluon 
propagator between quarks) always correspond to linear combinations of fully connected 
diagrams, (i.e. diagrams where all gluons are connected to the same quark-line) and are 
the only diagrams presently included in major event generators. When considering only 
physical topologies the ratio of fully connected to disconnected graphs changes slightly to 
the advantage of the fully connected graph, as the factor (l/2)# rin S s buildin S U P the tensor 
hits the disconnected topologies harder. 

The norm of the color tensor with all gluons attached to a quark going around in one 
direction is given by 



2N C 



J [(iV c 2 -l)^ + (-l)^(7V c 2 -l) 



The physical tensors, being sums of gluons attached to rings with quarks going around in 
opposite directions, contain mixed terms as well, these are however relatively suppressed, 



and for large N g or large N c Eq. (2.8) is a good approximation. Note that the norm grows 
as N c 9 , which is to be expected considering the N c = oo limit. It turns out, however, that 
it is easier to stick to the non-normalized versions of the color tensors. 

2.4 The case of both quarks and gluons 

In the general case of both external quarks and gluons the basis may be constructed by: 
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(1) Connect the quark lines to each other in all possible ways, giving N q \ possibilities. 

(2) For i = 1, 2, ...N g — 2, N g attach i of the N g gluons to the quarks in all possible ways. 

(3) Connect the remaining N g — i gluons as in the gluons only case, but keep cyclic and 
anti-cyclic tensors separately. 

(4) Distribute the quark and gluon indices in all possible ways among the different group- 
ings. 

The number of color tensors in this case grows slower than (N g + N q )\ but faster than 
(N g + N q )\/e, again giving a factorial growth. 

3. Calculating the effect of gluon exchange 

Below, the computational rules for gluon exchange will be derived, and it will be seen that 
exchanging a gluon trivially gives an explicit linear combination of the basis tensors. There 
is thus no need to calculate scalar products of the resulting color structure after exchange, 
with the basis tensors. As calculating scalar products was the most cumbersome part in 
previous calculations, this represents a major improvement. 

3.1 Computational rules 

In this section the computational rules for gluon exchange between the basis tensors con- 
structed in section || are derived. 

Note that the quarks in closed quark loops are just products of the way of writing 
down the basis and not physical particles, a gluon is thus never exchanged between the 
quarks in closed quark loops. 

We also have to decide on a convention for the triple gluon vertex. The convention 
used is 

f ei g with 

e = the external (incoming or outgoing) eikonal gluon index 

i = the internal (incoming or outgoing) eikonal gluon index (3-1) 

g = the soft exchange gluon index. 

This convention has the advantage that the sign is independent of how the diagram is 
drawn on a paper and whether a parton is incoming or outgoing. 

3.1.1 Gluon exchange between two quarks or anti-quarks 

In the simplest case a gluon is exchanged between two external quarks q\ and which 
in general have n and m gluons attached respectively. Using Eq. (p.2j), the effect of gluon 
exchange between the quarks q\ and qi in two different open quark lines may be written 

( <hffii-5ln<?i ® \ _^ 1 f ?iflli— 0in?2 ®\ J_ ( Ql9n-9lnqi <8> 

I <Z2S21-52m<?2 / 2 \q 2 g2i -9% m q\ ) 2N C I (Z 2 521-52m<?2 
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where the notation 

(Ql9ll---9lnQl ® \ _ ,gn ,gi2 +921 +922 f g2m /o o\ 

Q2921---92m,q2 J ~ " 5l dl20!l1 " d21 «2 <*22<*21— W2m-1 

is used. 

If the gluon is instead exchanged between the external anti-quarks, the indices on the 
quarks are kept whereas the indices on the anti-quarks are exchanged. 

3.1.2 Gluon exchange between quark and ant i- quark 

Exchanging a gluon between a quark q\ and an anti-quark q 2 results in 



qi9n-9i n qi ® \ _^ 1 / qi9ii--9in92i--92mq2 ®\ \_ I qi9n-9i n qi 

q292i--92 m q2 I 2\g 2 gi / 2N C \ q 2 g%\...g% m q% 



(3.4) 



The case where the involved quark and anti-quark are part of the same quark line can be 
obtained by identifying q 1 and q2 above. 

3.1.3 Gluon exchange between quark and gluon 

To derive the effect of gluon exchange between a quark and a gluon we use the relation 



Eq. (2A) to rewrite the triple gluon vertex. After this Eq. (|2.2| ) is applied (and it is noted 
that the 1/N C suppressed terms drop out). For gluon exchange between the quark q\ and 
the gluon g2i the result is: 

/ q 1 g u ...gi n qi ® \ _^ _1 / q 1 gii...gi n g2i+i-~g2 m q2 <8> | 1 ( q\g\\--g\ng2i--g2 n q2 ®\ 
\q~2921-g2i-g2mq2 J 2 ^2521-- -52^1 J 2 \ q 2 g2ig2i-iq\ J 

(3.5) 

If, in the left hand side above, g^i is in a closed quark loop this is accounted for by identifying 
q 2 and q2 , and if the gluon g2% is attached to the same quark-line as q\ this is taken care of 
by identifying q l and q2- 

3.1.4 Gluon exchange between anti-quark and gluon 

Employing the same calculational method as for qg results in 



/ <7iS , ii-"S , in<?i ® \ 1 / q\g2i---g2 m q2 ® \ _ l / qig2i+i-~g2 m q2 ® \ 
\l292i---g2i--g2m(l2 ) 2 y q 2 g 2 \ ...g2i-\g\\ ■■■g2 n qi J 2 \q 2 g2i---g2igii---gi n qi ) 

(3.6) 

where again, if g2% initially is in a closed quark loop this is accounted for by identifying q 2 
and q2, and if the gluon g2i was originally placed on the same quark line as g 1 this is taken 
care of by identifying q 2 and q\ . 
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3.1.5 Gluon exchange between two external gluons 

To derive the effect on the basis vectors of exchanging a gluon between two external gluons, 
two triple gluon vertices have to be replaced using Eq. fl2.1| ) and three gluon propagators 
have to be contracted using Eq. fl2,2|) . Again the non- leading N c terms drop out and the 
result of exchanging a gluon between gu and g2j is 



q\9x\-9u-9\nq\ ® \ _^ 

12921— 92j-g2mQ2 J 



(3.7) 



_i / qi9ii---gii-ig2j+i---g2 m q2 <s> \ 1 / q\9\\-9\i-\9\i92o+\-92mq2 ® j 
2 \~q~2g2i---g2j-ig2jgiigii+i---ginqi J 2 y<z 2 92i-"52j-i52j5ii+i"-5in<?i / 

+ l [ qi9ii-gu-i92j92j+x-92mq2 ® j _ l i qigii--gii-igug2j92j+i-92mq2 <s>\ 

2 \?2521-52i-l5li5li+l-5ln?l / 2 \ ^2521 ■■■92j-igii+l •••5ln9l / 

If one (or both) quark lines is (are) closed, then the corresponding quarks are to be iden- 
tified. If both gluons are part of the same quark line, then identify gig 2 , and (?29i if the 
quark line is closed. 

4. Some explicit examples 

4-1 gg -> 55 

As an explicit example of how the above strategy simplifies the problem of keeping track 
of the color structure, the process of 5152 — * 5354 will be considered in detail here. The 
soft anomalous dimension matrix for this case was first calculated in [12] and later, more 
elegantly in [14]. 

4.1.1 Construction of the basis 



To construct the basis the recipe outlined in section 2.3 is followed, starting with finding 
all the ways of grouping the gluons: 

(1) The four gluons can be grouped two and two {2, 2} or all four together {4}. 

(2) When all four gluons are attached to the same quark line, {4}, the indices can be 
placed in (4 — 1)! = 6 different ways. However, due to the symmetry between quarks 
and anti-quarks, clockwise and anti-clockwise gluon rings only enter in one linear 
combination, giving three physical tensors: 



^1929394 = (91525354) + (54535251) = Tr[WW*] + Tr[W«9i] 
^1929394 = (91525453) + (53545251) = TV[WW 3 ] + Tr[^W92 t 9l] 

<l m> = (51535254) + (54525351) = Tv[t^H^]+Tv[t^H^}. 

(4.1) 
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(3a) For the grouping {2,2}, the index order in the subgrouping doesn't matter (since 
Trf^tg] = T r [^2^i])- Each subgrouping thus only gives rise to one physical tensor. 

(3b, c) The gluon indices 9i, 92, 93, 94 may be split into the subgroupings as {{91, 92}, {93, 94}}, 
{{91, 93}, {92, 94}} and {{91, 94}, {92, 93}}, giving three basis tensors 



x x2 



Cg^sgi = (9i92)(9394) = Tv[t^]Tr[t^] = 1-1 <W<W 
C:; :m = {9i93){g2 9 A) = Tr[t^}Tr[t^n = S gig3 S g294 
C l 92939i = (9m)(g2 S 3) = Tv[t^]Tt[t 9 H^] = <WW (4.2) 



4.1.2 Calculation of soft anomalous dimension matrix 

As previously noted, once the effect of gluon exchange between g\ and 52 is calculated the 
effect of gluon exchange between any other gluons may be deduced. There may thus at 
most be ^basis different situations to keep track of. However, this number will in general be 
further reduced due to the irrelevance of non-participating indices. For example the effect 
of gluon exchange between g\ and 52 on Cg ig2g . ig4 is the same as the effect on Cg l92g3S4 - 
Thus, the physically different situations are: 

(1) A gluon is exchanged between two neighboring gluons on a quark ring with four 
gluons attached (for example gluon 1 and 2 on C 1 ). 

Applying Eq. (|3.7[) to the first half of Ch „ „ „ with the identification g\ — > g\2 = gu, 



31929394 

92 -> 921 = g2j, 93 -> 922 , 94 -> 9n and q 1 =q 2 ,q 2 = qi gives 



1 N c 
(91929394) -> -^(91 92) (9394) - -^(91929394)- (4.3) 



Similarly application to the second half results in 



1 N c 
(91949392) -> -^(91 92) (9394) - ^(91949392), (4.4) 



and it may be concluded that 



91929-3,94 2 91929394 2 9l929394 ' \^- u J 



(2) A gluon can be exchanged between two next to neighboring gluons. In this case we 
get for an exchange between g\ and 52 on C gig2g3gi 

C 3 — > C 5 -i- (d fi^l 

91929394 ^91929394^^91929394- v^- u ; 
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(3) A gluon may be exchanged between the gluons attached to a two gluon ring, such 
as g\ and g 2 on C gig2g . igi . This just gives a factor N c multiplying the old tensor, for 

y 9i929394 



example for gluon exchange between g\ and g 2 on C 4 



r* 4 — > N C A (A 7) 

31929394 c ^9l929394 - v '/ 



(4) A gluon may be exchange between two gluons attached to different two gluon rings 
such as gi and g 2 in C 6 gig2gagv giving 



r 6 

91929394 



-r 1 + -r 3 

2 91929394 ^ 2 91929394 " 



(4.8) 



The above information may be combined into a matrix describing the color algebra 
part for gluon exchange between g\ and 52 



pl2 

99^99 



I 2 




-1 


Vo 





-AT, 
2 














1 



1 
1 







2 

1 1 

2 2 

1 N r 








As the gluons 53 and 54 have the same relationship to each other in the basis as 
gi and g 2 , the color structure of the soft anomalous dimension matrix will be the same 



r<34 
^99 



^99- 



Similarly C 



->99 99^99 

The contributions C 14 



99^99 
and C^ 4 



r>23 



►99 alld C 99' 



^99^ 



99 99^99' 

may be calculated by using the results in Eqs. 
(4.5-|4.8) and relabeling indices. Letting T = O12 + ^34, U = ^13 + f^24 and V = ^14 + ^23 
be the phase space integrals the result can be written 



99^99 

1 _i 
2 



(4- 



V 



N C (T + V) 



U — T 


u-v 





■\N C (T + U) 

-\N C (U + V) 
V -T 
V-U T-U 
T — V 



u-v 
2 

2 










V-T 

2 

T-V 
2 



-AT C (7 




E/-T 
2 



T-U 
2 






To obtain physical results the scalar product matrix 
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99^99 — 

( AT c 6 -3AT c 4 +8Ar c 2 -6 



(4.10) 



-Nf+iN^-3 



-N*+4N*-3 j^-l) 2 l-N 2 c (iVg-l) 2 \ 



V 





4JV* 


±N'£ 


8N C 


8N C 


8N C 


-JV*+4iV"f-3 




-7V c 4 +47V 2 -3 






1-iVf 


4N£ 




iN'i 


8N C 


SN C 


87V C 


-JV*+4iVf-3 




N*-3N*+8N%-6 


1-JVf 




(Nl-lf 




W£ 


8N* 


8iV c 


8N C 


8N C 






1-iV 2 


(^-i) 2 


{n*-i) 


{N?-l) 


8N C 


8N C 


8N C 


16 


16 


16 


1-JV 2 




{Nl-lf 




(iv 2 -i) 2 


(Jv 2 -i) 


87V C 


8N C 


8N C 


16 


16 


16 


(^ c 2 -l) 2 


1-iVf 




K 2 -i) 


(^v 2 -i) 


(iv 2 -i) 2 


8JV C 


8N C 


8N C 


16 


16 


16 



calculated using Eq. (|Q|), is also needed. This matrix contains iV^^g 2 entries, however, 
closer consideration reveals that only six of them correspond to different contractions. 

It is worth remarking on the leading N c behavior of Eq. (fl.Sj), The computational 
rules in Eqs. (^2-3.7) contain no positive power of N c . Thus the N c in Eq. ( |4.9[ ) must 
come from closed quark loops. The only way to get a closed quark loop is to exchange a 
gluon between two neighboring partons attached to the same quark-line, i.e. only "color 
neighbors" radiate in the N c — > oo limit. The result after exchange contains a factor N c 
multiplying the old color structure. Leading N c contributions will therefore always be 
diagonal in the present basis. 

This is in close resemblance with the Dipole Cascade Model and the original Ariadne 
program in which only neighboring pairs of partons, dipoles, radiate [47-50]. In particular 
for gluon radiation from e + e~ — ► qq the leading N c piece should come from neighboring 
partons. It is cautioned, however, that there are many more non-leading N c contributions, 
than leading pieces, as there (in general) are many more non-neighboring partons. 

It is also worth remarking that although there are scalar products between different 
basis tensors that are suppressed by only one power of iV c , these scalar products are never 
between different fully connected topologies, i.e. never between tree level QCD gluon 
amplitudes. 



4.2 qq 



ma 



As an example of a process containing both quarks and gluons we consider the color struc- 
ture needed for gluon resummation for q\q 2 —>■ ftftgs- This color structure is important 
for (among other things) QCD corrections to the production of W^s decaying leptonically 
and being accompanied by three jets [51]. 



The result of constructing color tensors as outlined in section 2.4 is 



C liq2q 3 qA96 = (?4fl5«0 (ft ft) = *g?4<W 

C qi mi 95 = («49S) (ft 5592) = <W^291 

C iq2 q sq4gs = («4<75ft) (ft ft) = *gg 4 <W 

C qxqiqzq^ = feftXftSSft) = <W*a ft - ( 411 ) 
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Exchanging gluons between the partons in all possible ways results in a leading color 
diagonal part, al/N c suppressed off-diagonal part and a 1/-/V 2 suppressed diagonal part: 



A qq^qqg ~ 
N c 

— -Diagonal[fii 2 + fi 35 



( 



O15 



n 23 

^24 



45, 



O25 
"25 



^15 + ^25 + ^34, O 



'13 



25 



45, 



-flis + tt 2A + ^35] 



r>13 
^23 



Q 12 - 
Q u 

"14 — ^35 — ^45 
f^24 + ^35 + ^45 



"15 + ^23 
"15 + ^34 






^35 
^35 



^12 
^23 



S2i4 — O25 
il25 + ^34 







1 

"2W, 



n 



13 



O24 + 34 ) Diagonal[l, 1, 1, 1]. 



(4.12) 



Again, as we are working in a non-orthogonal basis, all scalar products are needed 



Jqq^qqg 



( \N C (iV 6 2 - 1) 

\N C (iV c 2 - 1) 

1) \ W " 1) 



i(iV 2 -l) i(^-l) \ 



3 - 1) 

IiV c (iV 2 - 1) 



- ,iV 2 

2 V c 







(4.13) 



iiV c (iV c 2 -l)/ 



4.3 grg -> 555 

The other color structure relevant for W plus three jets is that of qq — ► ggg. In this case 
an 11-dimensional matrix is needed to describe the color space (reducing to 10 for N c = 3). 
These results have been calculated and are electronically attached to this submission. Again 
there is a diagonal leading N c part, an off-diagonal part with relative suppression 1/N C and 
a 1/iV 2 suppressed diagonal contribution. 

5. Conclusions 



Q45 \ 
^45 



In the present paper a general recipe for constructing bases capable of dealing with the color 
structure needed for resummation for any number of colored partons has been presented. 
This in itself is a step forward. In addition the suggested bases are argued to have relatively 
nice computational properties. The bases are obtained from the N c = 00 case by splitting 
gluons in qq pairs and connecting color lines in all possible ways. The bases thus constructed 
will therefore neither be normalized or orthogonal for N c = 3, but they will span the 
space and have the property that gluon exchanges between any pair of external partons 
directly, i.e. without taking scalar products, result in linear combinations of basis vectors. 
Furthermore, as can be seen from the computational rules in Eqs. ( |3.2[]3.7| ), the result after 
gluon exchange contains at most four (often two or one) basis vectors, giving relatively 
sparse soft anomalous dimension matrices. 

The fact that there is no need to calculate scalar products to decompose the tensors 
resulting after gluon exchange is a major advantage. Otherwise there would, for each of the 
Np(Np - l)/2 possible gluon exchanges, be ~ A^asis ~ (^P-) 2 sect i° n liMH-l) sca l ar 
products to calculate. In addition, the computational time for calculating the effect of gluon 
exchange is further reduced, as the constructed bases maximally exploit the symmetry w.r.t. 
external parton indices. All indices corresponding to the same kind of parton enter the 
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basis on equal footing. Therefore, for example, once the effect of gluon exchange between 
any pair of gluons has been calculated, the effect of gluon exchange between any other can 
be obtained by relabeling of indices, corresponding to a renumbering of tensors. One thus 
at most has to calculate six (gg, qq, qq, or qq to qg and qg) different contributions to Y. 
For a hard scattering amplitude with only gluons it is enough to calculate one. 

From this it is clear that the major computational effort lies in computing the ~ -^basis 
scalar products between all the basis vectors. So far, it thus looks as if the calculational 
effort is reduced by a factor ~ iVp(iVp — l)/2 (times a factor coming from the fact that 
the basis vectors are simpler to take scalar products of). However, also in the case of 
calculating scalar products between basis tensors, the equal footing of the indices comes to 
rescue. This is so, as only the topology of the contraction, and not the labeling of indices, 
is important for determining the scalar product. 

Apart from the nicer scaling properties, the suggested bases have the advantage of 
being orthogonal and giving rise to diagonal soft anomalous dimension matrices in the 
N c — > oo limit. This enables a more straight forward comparison to event generators, 
tending to keep only the leading N c contribution. 
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